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Abstract 

We investigate the A(p, K + )^B reaction within a covariant two nucleon model. We 
focus on amplitudes which are described by creation, propagation and decay into 
relevant channel of A r *(1650), A r *(1710), and iV*(1720) intermediate baryonic res- 
onance states in the initial collision of the projectile nucleon with one of its target 
counterparts This collision is modeled by the exchange of tt as well as p and u 
mesons. The bound state nucleon and hyperon wave functions are obtained by solv- 
ing the Dirac equation with appropriate scalar and vector potentials. Expressions 
for the reaction amplitudes are derived taking continuum particle wave functions 
in both distorted wave and plane wave approximations. Detailed numerical results 
are presented in the plane wave approximation; estimates of the effects of the initial 
and final state interactions are given in the eikonal approximation. Cross sections 
of 1 - 2 nb/sr at the peak positions, are obtained for favored transitions in case of 
heavier targets. 

PACS numbers: 25.40.Fe, 13.75. - n„ 13.75. Jz 

Key words: Strangeness production, proton nucleus collisions, covariant two 
nucleon model. 



1 Introduction 



A hypernuclei are the most familiar hypernuclear systems which have been 
studied extensively by stopped as well as in-flight (K~ , 7r~) reaction [1,2,3] and 
also by K + ) reaction [4,5]. The kinematical properties of the (K~,tt~) 
reaction allow only a small momentum transfer to the nucleus (at forward 
angles), thus there is a large probability of populating A-substitutional states 
(A assume the same orbital angular momentum as that of the neutron being 
replaced by it). On the other hand, in the (tt + ,K + ) reaction the momentum 
transfer is larger than the nuclear Fermi momentum. Therefore, this reaction 
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can populate states with the configuration of an outer neutron hole and a A 
hyperon in a series of orbits covering all the bound states. During past years 
data on the hypernuclear spectroscopy have been used extensively to extract 
information about the hyperon-nucleon interaction (which could be quite dif- 
ferent from the nucleon-nucleon interaction) within a variety of theoretical 
approaches (see, e.g., [6,7]). 

Alternatively, A-hypernuclei can also be produced with high intensity pro- 
ton beams via p + A(N, Z) -> A B(N -l,Z)+n + K+, p + A(N, Z) -> 
A B(N,Z - 1) + p' + K + , and p + A(N,Z) -> A B(N,Z) + K+ reactions 
where N and Z are the neutron and proton numbers, respectively, in the 
target nucleus. In this paper we study the last reaction [to be referred to as 
A(p, K + ) A B] where the hypernucleus A B has the same neutron and proton 
numbers as the target nucleus A, with one hyperon added. This reaction is 
exclusive in the sense that the final channel is a two body system. First set 
of very preliminary data have already been reported for this reaction on deu- 
terium and helium targets [8]. More measurements for this reaction involving 
also the heavier targets are expected to be performed at the COSY facility of 
the Forschungszentrum Jiilich [see, e.g., Ref. [9]]. This reaction involves much 
larger momentum transfer to the nucleus as compared to the K + ) reac- 
tion [more than 1.0 GeV/c (see Fig. 1) as compared to only 0.33 GeV/c at the 
outgoing K + angle of 0°]. Therefore, it samples bound state wave functions 
in a region where they are very small and are unlikely to be reached in other 
reactions. From the spectroscopic point of view, the states of the hypernucleus 
A B excited in the (p, K + ) reaction may have a different type of configuration 
as compared to those reached in the K + ) reaction. Thus a comparison of 
informations extracted from the study of two reactions is likely to provide a 
better understanding of the hypernuclear structure. 

Theoretical studies of the A(p, K + ) A B reaction are rather sparse and pre- 
liminary in nature [10,11,12,13]. They are based on two main approaches; 
the one-nucleon model (ONM) [Fig. 2(a)] and the two-nucleon model (TNM) 
[Figs. 2(b) and 2(c)]. In the ONM the incident proton first scatters from the 
target nucleus and emits a (off-shell) kaon and a A hyperon. Subsequently, the 
kaon rescatters into its mass shell while the hyperon gets captured into one of 
the (target) nuclear orbits. Thus there is only a single active nucleon (impulse 
approximation) which carries the entire momentum transfer to the target nu- 
cleus. This makes this model extremely sensitive to details of the bound state 
wave function at very large momenta where its magnitude is very small lead- 
ing to quite low cross sections. In the ONM calculations of (p, K + ) and also of 
(p, 7r) reactions the distortion effects in the incident and the outgoing channels 
have been found to be quite important [13,14,15]. 

In the two-nucleon mechanism, on the other hand, the kaon production pro- 
ceeds via a collision of the projectile nucleon with one of its target counter- 
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Fig. 1. Center-of-mass momentum transfer available to the (p, K + ) reaction as a 
function of the laboratory bombarding energy (E p ) for 12 C and 40 Ca targets. For a 
fixed energy the allowed kinematic region extends from 0° to 180°. We have shown 
here only a few angles. 



parts. This excites intermediate baryonic resonance states which decay into a 
kaon and a A hyperon. The nucleon and the hyperon are captured into the 
respective nuclear orbits while the kaon rescatters into its mass shell. In this 
picture there are altogether three active bound state baryon wave functions 
taking part in the reaction process, thus allowing the large momentum transfer 
to be shared among three baryons. Consequently, the sensitivity of the model 
is shifted from high momentum parts of the bound state wave functions (not 
very well known) to those at relatively lower momenta where they are rather 
well known from (e, e'p) and (7,p) experiments and are relatively larger [see, 
e.g., [16]]. This could lead to larger cross sections. Moreover, in the TNM stud- 
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Fig. 2. Graphical representation of one-nucleon (a) and two-nucleon [(b) and (c)] 
models. The elliptic shaded area represent the optical model interactions in the 
incoming and outgoing channels. 

ies of other large momentum transfer reactions like (p, tt), (7, 7r), and (7,77) 
the distortion effects have been found to be less pronounced [17,18]. 

Most of the previous calculations [10,11,13] of the A(p, K + ) A B reaction have 
been done within the non-relativistic framework. However, for processes in- 
volving momentum transfers of typically 400 MeV/ c or more, a non-relativistic 
treatment of corresponding wave functions is questionable as in this region the 
lower component of the Dirac spinor is no longer negligible in comparison to 
its upper component (see, e.g, Ref. [19]). Furthermore, in the non-relativistic 
description one must resolve the ambiguities of the pion-nucleon-nucleon ver- 
tex that result from the non-relativistic reduction of the full covariant pion- 
nucleon-nucleon vertex (see, e.g., [17,19,20] for more details). This is impor- 
tant as this vertex is used to describe the initial collision between the incoming 
proton and the target nucleus. 

In this paper, we study the A(p, K + )\B reaction within a fully covariant TNM 
by retaining the field theoretical structure of the interaction vertices and by 
treating the baryons as Dirac particles. The initial interaction between the 
incoming proton and a bound nucleon of the target is described by the ir, p and 
uj exchange diagrams which were also used to describe the K + production in 
elementary proton-proton (pp) collisions [21,22,23]. In a previous study of the 
(p, K + ) reaction on a 40 Ca target [25] the latter two meson exchanges were not 
included. Some authors, however, have used an alternative approach for the 
elementary pp — > pK + A reaction [26,27,28] in which kaon exchange mechanism 
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between the two initial state nucleons leads to the K + production. We have not 
considered this mechanism here. There are indications, from from some results 
reported by COSY-TOF and ANKE groups [29,30], that the Dalitz plots for 
the pp — > pK + A reaction are dominated by the excitation of nucleon isobars, 
though modified by the Kp final state interaction (FSI) (see also Ref. [31]). 

In our model, the initial state interaction of the incoming proton with a bound 
nucleon leads to 7V*(1650)[§~], iV*(1710)[| + ], and iV*(1720)[§ + ] baryonic res- 
onance intermediate states which make predominant contributions to elemen- 
tary pp — > pK + A cross sections in the beam energy regime of near threshold to 
10 GeV [21]. Terms corresponding to the interference among various resonance 
excitations are included in the total reaction amplitude. We have ignored the 
diagrams of type 2(a) since contributions of such processes are expected to be 
very small in comparison to those of the TNM. 

In section II, we present the details of our formalism for calculating amplitudes 
corresponding to the diagrams shown in Figs. 2b and 2c, using continuum wave 
functions within both distorted wave (DW) and plane wave (PW) approx- 
imations. In section III, numerical results are presented for 4 He(j9, K + )\He, 
12 C(p, K+yftC, and 40 Ca(p, K + ) 4 iCa reactions within the the PW approxi- 
mation. We also give here the estimates of the effects of the initial and fi- 
nal state interactions which are calculated within the eikonal approximation. 
Summary, conclusions and future outlook of our work are given in section IV. 
Finally some formulas for the amplitudes are summarized in appendix A. 



2 Covariant two-nucleon model for A(p, K + ) reactions 



2. 1 Kinematics, elementary processes and contributing diagrams 



The elementary associate production reaction pp — * pK + A has been studied 
extensively both experimentally and theoretically (see, e.g., Ref. [32] for a 
recent review). The laboratory threshold energy for this reaction is 1.582 GeV. 
For the case of the A(p, K + )^B reaction the threshold energy is given by 



vthres ( m B- m A) , ( m K ~ ™n) , rn B m K 

tip = ~ 1 ~ 1 m N , [L) 

F IrriA 2m A txia 

where m B and are the masses of the hypernucleus \B and the target 
nucleus A, respectively and is the mass of the kaon. E l p hres is 0.739 GeV 
and 0.602 GeV for 12 C and 40 Ca targets, respectively which are considerably 
lower than the threshold energy value of the elementary production process 
as stated above. 
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The (p, K + ) reaction is characterized by the momentum transfer p cm = Pi—px 
to the residual nucleus. In Fig. 1, p cm are plotted as a function of the laboratory 
beam energy for several values of the outgoing kaon angles for 12 C and 40 Ca 
targets. It is clear that p cm varies between 1 GeV to 5 GeV. In the ONM 
a single nucleon alone carries this huge momentum to the nucleus. In the 
PW approximation, the ONM cross section of the (p, K + ) reaction is directly 
proportional to the modulus square of the bound state wave functions in the 
momentum space. Since at these large momenta the magnitudes of the bound 
state wave functions are very small, the corresponding cross sections are quite 
weak. Of course, the inclusion of the distortions changes this model towards 
a multi-nucleon mechanism thus shifting the sensitivity of the model to lower 
momentum components of the bound state wave functions. Therefore, apriori, 
these effects are very important in the ONM description of such reactions. 

The data on the associate production reaction pp — > pAK + , are available for 
beam energies ranging from very close to the threshold to upto 10 GeV which 
can be well described within an effective Lagrangian approach [21,22,23]. In 
this model, the associated kaon production proceeds via excitation, propaga- 
tion and decay of the iV*(1650), iV*(1710) and iV*(1720) intermediate baryonic 
resonant states, in the initial collision of two nucleons in the incident chan- 
nel. The coupling constants at the nucleon-nucleon-meson vertices are fixed 
by fitting to an elastic NN T matrix while those at resonance-nucleon-meson 
vertices are determined from the experimental branching ratios for the decay 
of the resonance into relevant channels except for those involving the u meson 
where they are determined from the vector meson dominance (VMD) hypoth- 
esis. To describe the near threshold data the FSI effects in the final channel 
are included within the framework of the Watson-Migdal theory. As an exam- 
ple, we show in Fig. 3 a comparison of calculated and experimental total cross 
sections for the pp — > pK + A reaction. The FSI amplitude has been calculated 
using the scattering length and the effective range parameters of the model A 
of A — p interaction described in Ref . [24] . The cross sections are plotted as a 
function of the excess energy (defined as ^/s — — — m\, where s/s 
is the invariant mass). It is clear that the model provides a good description 
of the experimental data in the entire range of beam energies. While at near 
threshold beam energies the reaction proceeds predominantly via excitation 
of the iV*(1650) resonance, the other resonances become more and more im- 
portant with increasing beam energies. One pion exchange graphs dominate 
the production process for all the energies. 

The structure of our TNM for the (p, K + ) reaction is similar to that of the ef- 
fective Lagrangian approach described above. The initial interaction between 
the incoming proton and a bound nucleon of the target is described not only 
by the dominant one-pion exchange mechanism but also by the p and oo ex- 
change processes. The latter may become more important at backward angles 
of the outgoing K + spectrum where momentum transfers are relatively larger. 
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Fig. 3. Comparison of the calculated total cross sections for the p + p — ► p + K + + A 
reaction with the corresponding experimental data [taken from Refs. [33,34]] as a 
function of the excess energy. The calculated cross sections represent the coherent 
sum of the amplitudes corresponding to all the contributing baryonic resonant states 
and the meson exchange processes. 



We use the same effective Lagrangians and vertex parameters to model these 
interactions. The initial state interaction between the two nucleons leads to 
the JV*(1650)[±~], 7V*(1710)[§ + ], and iV*(1720)[§ + ] baryonic resonance inter- 
mediate states. The vertex parameters here too are the same as those used in 
the description of the elementary reaction. Terms corresponding to the inter- 
ference between various amplitudes are retained. 

We denote the diagrams (and the corresponding amplitudes) by "projectile 
emission (PE)" if the intermediate meson starts from the projectile and rescat- 
ters from the target nucleon which is excited to an intermediate baryonic res- 
onance state [Fig. 2(c)] and by "target emission" if the opposite is the case 
[Fig. 2(b)]. In the two amplitudes there is a tremendous variation in the en- 
ergy of the intermediate meson. For example, in Fig. 2(b), the energy of the 
exchanged meson is given by the difference between single-particle binding 
energies of the initial and final intermediate states which is very small and 
can be approximated to zero. On the other hand, in Fig. 2(c), the energy of 
the intermediate meson is almost equal to the incident proton energy. Due 
to different four momentum transfers the behavior of the two amplitudes is 
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quite different and they explore different regions of the intermediate meson 
propagator. Thus corresponding medium effects have to be known over a wide 
range of energy and momenta in order to incorporate their effects into the 
calculations. 



2.2 Calculations of the TNM amplitudes 



The effective Lagrangian for the nucleon-nucleon-meson vertices are given by 



Xnn« = -f^757,r • (2) 

&NN P = -9NN P ^N (l» + T • P ^ N . (3) 

&NN<o = -9NN U *N (l» + ^^ a ^ d ^j (4) 

We have used notations and conventions of Bjorken and Drell [35]. In Eq. (2) 
Til jv denotes the nucleon mass. In Eqs. (3) - (4) is defined as 



<V = 2 (7^-7*7/,) (5) 

The A/W-meson vertices are corrected for off-shell effects by multiplying the 
corresponding coupling constants by form factors of the forms [21] 



J?NN _ ( m i \ f R \ 

where qi and are the four momentum and mass of the ith exchanged meson 
and \i are the corresponding cutoff parameters [their values are taken to be 
the same as those given in Ref. [21]]. The latter governs the range of the 
suppression of the contributions of high momenta. 

As the A-hyperon has zero isospin, only isospin-1/2 nucleon resonances are 
allowed. Below 2 GeV center of mass (cm.) energy only three resonances, 
AT* (1650), AT* (1710), and A 7 "* (1720), have significant decay branching ratios 
into the K + A channel. As in the description of the associate K + production 
in the elementary pp collisions, only these three resonances have been consid- 
ered in this work. The effective Lagrangians for the resonance-nucleon-meson 
vertices are written as [21,36,37,38,39] 

■&N; /a Nn = -9N* /2 N^N* 1/2 i^r^ + h.c, (7) 
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^n* 1/2 n p = -9n* 1/2 n p ^n* ^ V^STr ■ ff^ N . + h.c. (8) 

^N* /2 Nu> = ~9Nl /2 N^N* ^-T^d U UJ^ N . + h.C. (9) 

QN* Ntt - 

J?n; /2 n, = — — — ^^*FnT ■ + h.c. (10) 

6 J. IT) " 



9N* Np 

-N'Np = i : *pT> 

3 / 2 m N * + m N 



&n-n p = i 3J t V (0"p" - 0"p") 7,75^iv + h.c. (11) 



^• /2 iv w = i 9N;, f" (Vu» - d^u") 7,75^tv + h.c. (12) 

3 / 2 m N * + 

The operator r(r 7r ) is either 75 (unity) or unity (75) depending on whether 
the parity of the resonance is even or odd, respectively. The operator is 
Cftu or 7507^ for these two cases. ^ is the A r *(1720) vector spinor. It may be 
noted that we have used the conventional Rarita-Swinger coupling for the spin- 
| particle [36,37] which also has off-shell contributions to the spin-| partial 
waves. To describe the latter, an operator Q afl (z) = g au — |(1 + 2,2)7^7^ has 
also been included in the vector spinor vertex [36,38,39,40]. The choice of 
the off-shell parameter z is arbitrary; it is treated as a free parameter to be 
determined by fitting to the data. In Ref. [41] an alternative interaction has 
been suggested that does not activate the spurious spin-i degree of freedom. 
We, however, work with the Lagrangians as given in Eqs. (10)-(12). It may 
further be noted that we have used a pseudovector (PV) coupling for the NNk 
vertex and a pseudoscalar (PS) one for Ny 2 Nn and Ny 2 AK vertices. This 
provides the best description of the pp — > pAK + data [21]. 

The effective Lagrangians for the resonance-hyperon-kaon vertices are written 

as 



^N* /2 AK+ = -gNZ /2 AK+*N'irT& K +V + h.C. (13) 

9N* AK+ - 

&N-AK+ = — P^r • dP* K+ V + h.c. (14) 

3 / 2 1TIk+ 

Signs and values of various coupling constants have been taken from [21,25] 
and are shown in Table 1. It should be noted that for the nucleon-nucleon- 
meson vertices, the coupling constants were taken to be energy dependent in 
the same way as that described in Ref. [21]. The iVjy 2 iV-meson and N£, 2 N- 
meson vertices are corrected for the off-shell effects by multiplying the corre- 
sponding coupling constants by form factors of the forms (which is also used 
in Ref. [42]). 



NN* 



(\N*\4 



(Af )4 +(g 2_ m 2 )2 
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Table 1 

Coupling constants for various vertices used in the calculations. 



Vertex 


Coupling Constant (g) 


NNn 


12.56 


NNp 


2.00 


NNlo 


24.05 


N*(1710)Ntt 


1.04 


N*(l710)Np 


4.14 


N*(1710)Nu 


1.22 


N*(17W)AK+ 


6.12 


N*(W50)Nir 


0.81 


N*(W50)Np 


2.62 


N*(1650)Nu> 


1.80 


N*(1650)AK+ 


0.76 


N*(1720)Ntt 


0.21 


N*(l720)Np 


33.75 


N*(1720)Nu 


16.93 


N*(1720)AK+ 


0.87 



The values of the cut-off parameters appearing therein are taken to be 1.2 
GeV in all the cases in order to reduce the number of free parameters. It may 
noted that same form of F^ N * and the value of the cut-off parameter was also 
used in the calculations of the elementary cross sections shown in Fig. 3. 

After having established the effective Lagrangians and the coupling constants, 
one can write down, by following the well known Feynman rules, the ampli- 
tudes for graphs 2(b) and 2(c). The isospin part is treated separately which 
gives rise to a constant factor for each graph. For example, the amplitude for 
graph 2(b) with one-pion exchange mechanism and excitation of the positive 
parity spin-| baryonic resonance is given by, 



M 2b (N* 1/2 ) = Co(g^) (9N l/2 N„)(g N1/2 AK + MP2)w»q' 
x ip(p 1 )D 7T (q)'ijj(p A )-f 5 D N * (pjv*)75 
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(+)/ 



(16) 



where various momenta are as defined in Fig. 2(b). In addition, p^* is the mo- 
mentum associated with the intermediate resonance and p/^ is that associated 
with the A hyperon. The isospin factor Cf s 6 is unity. For the graph 2(c) also it 
is unity. The functions ip are the four component (spin space) Dirac spinor in 
momentum space [19,43]. {p'kiPk) ^^{p'iiPi)] is the wave function for 
the outgoing kaon [incoming proton] with appropriate boundary conditions. 
Function -Dj(pj) is the propagator for the ith particle with momentum pi. The 
propagators for mesons and spin-| and spin-| intermediate resonances are 
given by 




D N *(p)=i 



q 2 -m 2 -U p (q)) ' 

P v l v + m N * /2 

P 2 -(m N « /2 - i r N *j2y 

i($ + m N * /2 ) 
P 2 -(m N;/2 -iT N;/ j2) 2 



(17) 
(18) 

(19) 



3 3m 2 M * 

JV 3/2 



3 7YI AT* 

iV 3/2 



(//-//)]■ (20) 



In Eq. (17) II 7r (g) is the (complex) pion self energy which accounts for the 
medium effects on the propagation of the pion in the nucleus. Similarly, U p (q), 
in Eq. (18) is the same for the p meson. We have adopted the approach of 
Ref. [44] for the calculation of II 7r (g) which has been renormalized by including 
the short-range repulsion effects through the constant Landau-Migdal param- 
eter g'. All the relevant formulas for this calculation are given in Ref. [19]. Like 
(p, it) [19] and (p, p'ir) [45] reactions, (p, K + ) cross sections too are sensitive 
to the parameter g' [25]. For the self energies of p and u mesons, we have 
followed the procedure described in Ref. [46]. 

In Eqs. (19) and (20), Tn* is the total width of the resonance which is in- 
troduced in the denominator term to account for the finite life time of the 
resonances for decays into various channels. T N * is a function of the center of 
mass momentum of the decay channel, and it is taken to be the sum of the 
widths for pion and rho decay (the other decay channels are considered only 
implicitly by adding their branching ratios to that of the pion channel) [21]. 
The medium corrections on the intermediate resonance widths have not been 
included. We do not expect any major change in our results due to these ef- 
fects. As is pointed out in [47,48,49], the medium correction effects on widths 
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of the s- and p-w&ve resonances, which make the dominant contribution to the 
cross sections being investigated here, are not substantial. On the other hand, 
any medium modification in the width of the <i-wave resonance is unlikely to 
alter our results as their contributions to the cross sections are negligible [46]. 

The four component Dirac spinors ip{pj) are the solutions of the Dirac equation 
in momentum space for a bound state problem in the presence of an external 
potential field [19,50] 



fyip) = m N ip(p) + F(p), (21) 
where 



F(p)=6(po-E) 



|rfVK(- P ')^(p + p') 



7o / dWf-pWp + p') 



(22) 



In our notation 'p' represents a four momentum, and 'p' a three momentum. 
The magnitude of p (\p\) is represented by k, and its directions by p. 'po' rep- 
resents the time-like component of the momentum. Similarly, the magnitude 
and directions of the radial vector r is represented by a and f, respectively. 
In Eq. (22), V s and V® represent a scalar potential and time-like component 
of a vector potential in the momentum space. The spinors ip(p) and F(p) are 
written as 



ip(p) = 6(po-E)< 



Ftp) = 5tP«-E){ ™*™<*}\ (23) 

where f(k)[((k)] is the radial part of the upper component of the spinor 
ip(p)[F(p)]. Similarly g(k)[('(k)] are the same of their lower component. Func- 
tions f(p) and g{p) are the Fourier transforms of the radial parts of the cor- 
responding coordinate space spinors. C(^) s are related to /, g and the scalar 
and vector potentials in the following way 



C(*) = C<(*) -C«(*), 

C(k) = C(k)-C(k), (24) 



where 



12 



Uk) = 4,Jdk'k'W s (k')0±^- v (25) 

Uk) = 4vr / dk'k'^(k') ^ f ^ r (26) 

C(k) = Anjdk'k-V s (k')^±^, (27) 

= 4tt / dk'k'^k') ^^, . (28) 

More details of the derivations Eqs. (21)-(28) can be found in appendix A of 
the Ref. [19]. In Eq. (23), ^y^j are the coupled spherical harmonics 



&nk = < trntlfif* \jrrij > Y lmi (p)xi/2 W , (29) 

where £' = 21 — j with £ and j being the orbital and total angular momenta, 
respectively and Y represents the spherical harmonics. Xi/2 W is the spin space 
wave function of a spin-| particle. In this paper we consider only those cases 
where the final state has a pure single particle configuration. Likewise, the 
intermediate states are supposed to have pure single-hole or particle structures. 

The incident proton spinor and the outgoing kaon field are given by 

J p Lp Alp 

® { k ] *(p'k,Pk) = 5(p' K0 - E K ) £ (-Y«Y eKmK (p K )Y; KmK ip' K ) 

£ K m K 

xf EK (k' K ,k K ), (31) 

where Ei and E K represents the energies of the incident proton and out going 
kaon, respectively, and where we have defined 

oo 

Fl v j p {K k i) = ^J 3L p {k'fi) F th^ a ^ da > ( 32 ) 



1 00 

Gl> p j p (Ki ki) = ^2 j h' p {K°) G L' pJp (ki, a)a 2 da, (33) 
o 

1 °° 

fe K ( k ' K , k K ) = — J je K (k' K a) ff k {k K , a)a 2 da. (34) 
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In Eqs. (32) and (33), Ff 3 and G C L , j are the coordinate space spinors which 
are obtained by solving the Dirac equation for the scattering state with ap- 
propriate complex scalar and vector potentials. In Eqs. (30)-(35), the mo- 
mentum arguments p denote the quantum numbers of the asymptotic free 
states whereas p' represent the momentum coordinates. A short derivation of 
Eqs. (30)-(31) is presented in appendix B. In Eq. (34) wave function ff is the 
coordinate space solution of the Klein-Gordon equation with kaon-nucleus op- 
tical potential (see, e.g., Ref [51]). In these equations ji represent the spherical 
Bessel function. It should be mentioned here that due to oscillatory nature of 
the asymptotic these wave functions in the asymptotic region, the integrals 
involved in Eqs. (32)- (34) converge poorly. Such integrals can, however, be 
calculated very accurately by using a contour integration method [52]. 

In the PW approximation, the wave functions ^i + \p'i,Pi) and $k^*(p'k,Pk) 
are given by 



*? , V«*) = J B £r L ( *" )*«-*), (35) 
V 2m N V^+^W 

^k ) *(PK,PK)=S\p > K -p K ). (36) 

In Eq. (35) represents the two component Pauli spinor. The other symbols 
have the same meaning as in Ref. [35]. 

By repeated use of Eqs. (21) and (22) one can perform the tedious but straight 
forward algebra to trace out the 7 matrices from the expression for the ampli- 
tudes [e.g., Eq. (16)]. Furthermore, by using Eqs. (21)-(28) and (30)-(33) and 
performing the angular momentum algebra one can get the amplitudes in a 
form which is suitable for numerical evaluation. We write in appendix A, e.g., 
the final form of the amplitude M 2 b [Eq. (16)] in both full distorted wave as 
well as plane wave approximations. 

To get the T matrix of the (p, K + ) reaction, one has to integrate the amplitude 
corresponding to each graph over all the independent intermediate momenta 
subject to constraints imposed by the momentum conservation at each vertex. 
For instance, for the amplitude corresponding to Eq. (16) the respective T 
matrix is given by 



T (AT* \- f f dA PK f d A pA f d 4 p 2 

W^i,2)-j (2n)4 J (27T)4 J (2n)4 J {2 , 



(2tt) 4 7 (27r) 4 7 (2tt) 4 7 (2tt) 4 
xS(q-pi +P2)5(p N * -p K - Pa) 

x 6( PA -p' i -q + p' K )M 2b (N* /2 ). (37) 
It can be seen, from Eq. (37), that in the full distorted wave theory one would 
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be required to perform a twelve dimensional integration to evaluate the T- 
matrix T 2 } } which requires a tremendous numerical effort. This gets further 
complicated due to the fact that the number of partial waves required in 
Eqs. (30)-(32) are quite large for these high energy reactions. On the other 
hand, in the plane wave approximation the integrations over variables p\ and 
p' K become redundant. This not only reduces the dimensionality of the inte- 
grations by a factor of 2 but also removes the requirement of partial wave 
summations altogether. In this exploratory study we would like to make cal- 
culations for very many cases involving a variety of targets and beam energies 
in order to understand the basic mechanism of this reaction and to get the 
estimates of the order of magnitudes of various cross sections to guide the 
planning of the future experiments. We have, therefore, made our calcula- 
tions in a numerically simpler plane wave theory. However, estimates of the 
distortions have been made within the eikonal approximation. 

The differential cross section for the (p, K + ) reaction is given by 



dcr 1 m p m A m B p K ^ 2 

dn-(^r(E Pi + E A yp t m % llm ^ 1 ' [6 * } 

where E p . and E A are the total energies of incident proton and the target 
nucleus, respectively while m p , uia and m B are the masses of the proton, and 
the target and residual nuclei, respectively. The summation is carried out over 
initial (mj) and final {rrif) spin states. T is the final T matrix obtained by 
summing the transition matrices corresponding to all the graphs. 



3 Results and Discussions 



3.1 The nuclear and hypernuclear structure and bound state spinors 



The spinors for the final bound hypernuclear state (corresponding to mo- 
mentum p\) and for two intermediate nucleonic states (corresponding to mo- 
menta p\ and P2) are required to perform numerical calculations of various 
amplitudes. We assume these states to be of pure-single particle or single-hole 
configurations with the core remaining inert. To simplify the nuclear structure 
problem the quantum numbers of the two intermediate states are taken to be 
the same although it is straight forward to include also those cases where they 
may occupy different orbits. The spinors in the momentum space are obtained 
by Fourier transformation of the corresponding coordinate space spinors which 
are the solutions of the Dirac equation with potential fields consisting of an 
attractive scalar part (V s ) and a repulsive vector part (V v ) having a Woods- 
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q (fm" 1 ) 

Fig. 4. (Upper panel) Momentum space spinors [<&((/)] for Osi/ 2 A orbit in j(He 
hypernucleus. \F(q)\ and \G(q)\ are the upper and lower components of the spinor, 
respectively. (Lower panel) Momentum distribution [p(q) = \F(q)\ 2 + |G(g)| 2 ] for 
the same hyperon state calculated with these wave functions. 

Saxon form. This choice appears justified as the Dirac Hartree-Fock calcula- 
tions [53,54] suggest that these potentials tend to follow the nuclear shape. 
The same potential form has also been used in the relativistic ONM [15] and 
TNM calculations [19] of the (p, n) reaction. 

With a fixed set of the geometry parameters (reduced radii r s and r v and 
diffusenesses a s and a v ), the depths of the potentials were searched in order 
to reproduce the binding energies of the particular state (the corresponding 
values are given in Table 2). We use the same geometry for the scalar and 
vector potentials. The depths of potentials V s and V v were further constrained 
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Fig. 5. same as in Fig. 4 for 0p 3 / 2 A and 0si/ 2 A orbits in 



hypernucleus. 



by requiring that their ratios are equal to -0.81 as suggested in Ref. [55]. Ex- 
perimental inputs have been used for the configurations and the single baryon 
binding energies of the states in hypernuclei iHe [2] and [2,3]. However, 
the quantum numbers as well as the single baryon binding energies for states 
in 4 aC hypernucleus have been taken from the density dependent relativis- 
ts hadron field (DDRH) theory predictions of Ref. [7], which reproduce the 
corresponding experimental binding energies reasonably well. In this case, we 
have compared, for each state, the spinors calculated by our well-depth search 
method with those calculated within the DDRH [7] theory and find an excel- 
lent agreement between the two. 

In Figs. (4)- (6) we show the lower and upper components of the Dirac spinors 
in momentum space for the Osi/2 hyperon in jjHe, and Osi/2 hyperons 

in ^C, and 0g? 3 / 2 and Opi/2 hyperons in 4 lCa, respectively. In each case, we 
note that only for momenta < 1.5 fm -1 , is the lower component of the spinor 
substantially smaller than the upper component. In the region of momentum 
transfer pertinent to exclusive kaon production in proton-nucleus collisions, 
the lower components of the spinors are not negligible as compared to the 
upper component which clearly demonstrates that a fully relativistic approach 
is essential for an accurate description of this reaction. 

The spinors calculated in this way provide a good description of the experi- 
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Table 2 

Searched depths of vector and scalar potentials and the binding energies of the A 
and nucleon bound states. 



State 


Binding Energy (e) 


v v 


r v 


a v 


V a 


r s 


a s 




(MeV) 




(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


AHe(0ai /2 ) 


3.12 




135.424 


0.983 


0.578 


-167.190 


0.983 


0.606 


4 He(0s 1/2 ) 


19.814 




384.053 


0.983 


0.578 


-429.695 


0.983 


0.606 




0.860 




180.712 


0.983 


0.578 


-223.101 


0.983 


0.606 


1 iC(0p 1/2 ) 


0.708 




206.161 


0.983 


0.578 


-254.520 


0.983 


0.606 




11.690 




166.473 


0.983 


0.578 


-205.522 


0.983 


0.606 


12 C(0p 3 / 2 ) 


15.957 




382.598 


0.983 


0.578 


-472.343 


0.983 


0.606 


4 ACa(0s 1/2 ) 


17.882 




154.884 


0.987 


0.676 


-191.215 


0.982 


0.700 


4 ACa(0p 3 / 2 ) 


9.677 




179.485 


0.987 


0.676 


-221.587 


0.982 


0.700 


4 ACa(0 Pl/2 ) 


9.140 




188.188 


0.987 


0.676 


-232.331 


0.982 


0.700 


4 ACa(0d 5/2 ) 


1.544 




207.490 


0.987 


0.676 


-256.160 


0.982 


0.700 


4 ACa(ls 1/2 ) 


1.108 




192.044 


0.987 


0.676 


-237.091 


0.982 


0.700 


4 ACa(0(i3 /2 ) 


0.753 




230.206 


0.987 


0.676 


-284.205 


0.982 


0.700 


40 Ca(0d 3/2 ) 


8.333 




360.980 


0.987 


0.676 


-445.660 


0.982 


0.700 



mental nucleon momentum distributions for various nucleon orbits as is shown 
in Ref. [19]. In the lower panel of each of Figs. (4)-(6), we show momentum 
distribution p(q) [= F(q)\ 2 + |G(g)| 2 ] [16] of the corresponding A hyperon. It 
can be noted that in each case the momentum density of the hyperon shell, 
in the momentum region around 0.35 GeV/c, is at least 2-3 orders of magni- 
tude larger that around 1.0 GeV/c. Since, in the TNM the large momentum is 
shared among three baryons, the model is sensitive to the bound state spinors 
in the momentum regime where they are well described and are quite large. 
Therefore, one expects to get a larger cross section for exclusive meson pro- 
duction reactions within this model. 

3.2 Cross sections for A(p, K + )\B reaction 

3.2.1 Results within the plane wave approximation 

The self-energies U n , U p and Il w of pion, rho and omega mesons, respectively, 
are another input quantities required in the calculations of the kaon produc- 
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Fig. 6. Same as in Fig. 4 for 0d 3 / 2 A and 0pi/ 2 A orbits in 



4 ACa hypernucleus. 



tion amplitudes. They take into account the medium effects on the interme- 
diate meson propagation. The p and ui self energies have been calculated by 
following the procedure described in Ref. [46]. It may be noted that the inter- 
mediate mesons are always space- like and their self energies are functions of 
the exchanged four momenta. 

The pion self energy, which is more crucial as one-pion exchange diagrams 
dominate the (p, K + ) cross sections, is the same as that shown in Ref. [25]. 14^ 
is obtained by calculating the contribution of particle-hole (ph) and delta-hole 
(Ah) excitations produced by propagating pions [44] . This has been renormal- 
ized by including the short-range repulsion effects by introducing the constant 
Landau-Migdal parameter g' which is taken to be the same for ph — ph and 
Ah — ph and Ah — Ah correlations which is a common choice. The parameter 
g', acting in the spin-isospin channel, is supposed to mock up the compli- 
cated density dependent effective interaction between particles and holes in 
the nuclear medium. Most estimates give a value of g' between 0.5 - 0.7. The 
sensitivity of the (p, K + ) cross sections to the parameter g' is discussed in 
Ref [25] where a detailed discussion is presented of the effect of the self-energy 
correction (to the pion propagator) on various amplitudes. 

In the PE graph the values of the time-like momentum component go m the 
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e K (deg) 

Fig. 7. Relative contributions of target emission (TE) (solid line) and pro- 
jectile emission (PE) (dashed line) graphs to the differential cross section for 
the 40 Ca(p, K + ) 4 \C&(0d 3 / 2 ) reaction at the beam energy of 1.5 GeV. Only 
pion-exchange graphs are included in these calculations. The pion self-energy used 
in these calculations has been renormalized by using the Landau parameter g' = 0.5. 

meson propagator is taken to be go = Ei, where the latter is the energy 
of the incident proton (ignoring the nucleon binding energies) by the energy 
conservation at the vertices. Due to this the meson propagator in this diagrams 
[Fig. 2(c)] may develop a pole in the integration at \q\ = \J\q\ 2 + m\ {m\ is the 
mass of the meson). With the inclusion of the self-energy (which, in general, 
is complex at non-zero energies), the pole is automatically removed from the 
real axis and the pole integration problem disappears and the graph 2(c) 
becomes well behaved. Furthermore, this leads to a strong reduction in the 
contributions of this graph. 

On the other hand, such a pole does not develop in the TE diagram [Fig. 2(b)] 
since in this case the pion self-energy is real and attractive due to q ~ [as 
can be seen from Fig. 3 of Ref. [25]]. This leads to a small denominator of the 
pion propagator leading to a consequent increase in the contribution of this 
graph. In fact this amplitude dominates the (p, K + ) cross sections as can be 
seen from Fig. 7 where we show the relative contributions of graph 2(b) and 
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e K (deg) 

Fig. 8. Differential cross section for the 12 C(p,K + ) 1 ^C reaction for the incident 
proton energy of 1.8 GeV for various bound states of final hypernucleus as indicated 
in the figure. The A separation energies for various states were the same as those 
shown in Table 2. 

2(c) to this cross section considering only the one pion exchange mechanism. 
Similar domination of the TE amplitudes have also been seen in case of the 
TNM calculations of the (p, n) reaction [19]. 

In Figs. 8 and 9, we show the kaon angular distributions corresponding to 
various final hypernuclear states excited in reactions 12 C(p, K + ) l ^C and 
40 Ca(p, K + ) A \C&, respectively. The incident proton energies in the two cases 
are taken to be 1.8 GeV and 1.5 GeV, respectively where the angle integrated 
cross sections for the two reactions are maximum. In both the cases we have 
taken the Landau parameter g' = 0.5 in the calculations of the pion self en- 
ergy. The calculations are the coherent sum of all the amplitudes correspond- 
ing to the various meson exchange processes and intermediate resonant states. 
Clearly, the cross sections are quite selective about the excited hypernuclear 
state, being maximum for the state of largest orbital angular momentum. This 
is due to the large momentum transfer involved in this reaction. It may be 
noted that in each case the angular distribution has a maximum at angles 
larger than 0°. This is due to the fact that there are several maxima in the 
upper and lower components of the momentum space bound spinors in the 
region of large momentum transfers. Therefore, in the kaon angular distribu- 
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tion the first maximum may shift to larger angles reflecting the fact that the 
bound state wave functions show diffractive structure at higher momentum 
transfers. 

We see that while in Fig. (8) there is a difference of more than an order of 
magnitude between the cross sections for the p-shell and the s-shell excitations, 
the two are almost of the same order of magnitude in Fig. 9. To understand 
this, we note from Table 2 that for the hypernucleus, the binding energies 
of the p-shell states are in the range of only 0.70-0.90 MeV as compared 
to 11.69 MeV of the s-shell state. The bound state spinor (<3>) behaves, at 
larger momentum transfers, approximately as where k 2 oc e. Thus, with 
increasing binding energies the magnitudes of $ decreases at larger momenta. 
This leads to a decrease in the cross sections of those processes which are more 
sensitive to the larger momentum transfers. The p-shell state spinors in l \C, 
are considerably larger than those of the s state in the relevant momentum 
transfer (q) region (~ 2-4 fm _1 ). Therefore, thep-shell transitions are enhanced 
as compared to the s-shell one. Due to this binding energy selectivity, we call 
the reaction leading to the p-shell excitation as a matched transition while 
that to the s-shell a mismatched one. 

On the other hand, both p- and s-shell transitions in Fig. 9 are strongly 
mismatched due to the fact that for the 4 A.Ca hypernucleus both these states 
have large binding energies as compared to those of the ci-shell. In this case, the 
spinors of the p- and s-states do not differ substantially from each other in the 
relevant q region, which in turn makes the corresponding cross sections not to 
differ too much from each other. As a matter of fact, states with larger angular 
momenta are squeezed in the coordinate space (r-space) due to the centrifugal 
barrier and hence they are spread out in the momentum space (q-space) which 
enhances the cross sections for correponding transitions. It should, however, be 
mentioned here that for a fixed angular momentum, states with larger binding 
are more compact in r-space, and therefore are wider in q-space. Hence, in 
such cases cross sections to states with larger binding may actually be larger. 

The absolute magnitudes of the cross sections near the peak is around 1-2 
nb/sr, although the distortion effects could reduce these values as is shown 
below. This order of magnitude estimates should be useful in planning of the 
future experiments for this reaction. As found in Ref. [25] contributions from 
the iV*(1710) resonance dominate the total cross section in each case. We also 
note that the interference terms of the amplitudes corresponding to various 
resonances are not negligible. It should be emphasized that we have no freedom 
in choosing the relative signs of the interference terms. 

In Fig. (10), we show the kaon angular distributions for reaction 4 He(p, K + )\He 
leading to the 0si/ 2 state of the hypernucleus at the proton beam energy of 
2.0 GeV. In contrast to Figs 8 and 9, the maximum in the cross section, in this 
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Fig. 9. Differential cross section for the 40 Ca(p, K + ) A \C& reaction for the incident 
proton energy of 1.5 GeV for various bound states of final hypernucleus as indicated 
in the figure. The A separation energies for various states were as shown in Table 2. 

case, occurs at zero degree and it decreases gradually as the angle increases. 
This difference in the pattern of angular distribution seen in Fig. 10 from that 
of Figs. 8 and 9 can be understood from the fact that for momentum transfers 
relevant to this case the Dirac spinors are smoothly varying and are devoid 
of structures as can be seen in Fig. 4. In any case, from the purely kinemat- 
ical arguments it is clear that the maximum in the cross section for a £ = 
transition is expected normally to occur at smaller angles as compared that 
for a £ 7^ one. The absolute magnitude of the cross section at the forward 
angles, in this case, is about 0.02 nb/sr. This value is somewhat smaller than 
the upper limit of the experimental center of mass cross section deduced for 
this reaction in a very preliminary study made in Ref. [8] at a similar beam 
energy [8]. 

In Fig. (11)-(12), we investigate the role of various meson exchange processes in 
describing the differential cross sections as a function of outgoing kaon angles 
for reactions 12 C(jo, K+yiC(0p 3/2 ), and 40 Ca(j», K + ) A lCa(0d 3/2 ) at beam ener- 
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Fig. 10. Differential cross section for the 4 He(j>, i"C + )j{He reaction for the incident 
proton energy of 2.0 GeV for the bound state of final hypernucleus as indicated in 
the figure. The A separation energy was as shown in Table 2. 



gies of 1.8 GeV, and 1.5 GeV, respectively. In each case, dashed, dashed-dotted 
and dotted lines show the contributions of ir, p and uo exchange processes, re- 
spectively. Their coherent sum is depicted by the full line. We note that pion 
exchange graphs dominate the production process for all angles in each case. 
This result is similar to the observation made in the case of associated kaon 
production in elementary pp collisions [21,56]. For the 12 C(p, K + ) 1 ^C(0j?3/2), 
and 40 Ca(p, K + Y\Ca,{{)d^/2) reactions, contributions of p and u exchange pro- 
cesses are unimportant at forward angles. However, at backward angles they 
become relatively more significant. This is understandable as at the backward 
angles the momentum transfer is larger which leads to greater contributions 
from the heavier meson exchanges. We see a destructive interference between 
the 7r and p contributions. The p and uo exchange contributions put together 
do not affect appreciably the pion exchange only cross sections even at the 
backward angles due their mutual cancellation. 



In Figs. 13 and 14 we show the predicted beam energy dependence of the angle 
integrated cross sections for 12 C(jo, K + ) 1 \C(0p 3 /2) and 40 Ca(p, -fT + ) 4 ACa(0<i3/2) 
reactions, respectively. The results indicate a substantial dependence on the 
beam energy of the integrated cross sections. As the beam energy increases be- 
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Fig. 11. Differential cross section for the 12 C(p, K + ) 1 ^C(0p 3 / 2 ) reaction for the inci- 
dent proton energy of 1.8 GeV. The dashed, dashed-dotted and dotted line show the 
contributions of ir, p and u> exchange processes, respectively. Their coherent sum is 
depicted by the full line. 

yond the respective threshold energies, Ep hres , the cross section first increases 
fast and then after peaking around a given value starts decreasing slowly. This 
trend is reminiscent of the general nature of the beam energy dependence of 
associate kaon production cross sections in elementary pp collisions. The slight 
shift in the peak position towards the smaller value in the case of the heav- 
ier nucleus reflects the fact that the value of E l ^ res is smaller as one goes to 
heavier target. 

3. 2. 2 Effects of initial and final state interactions 

Calculations presented thus far have been done in a plane wave approximation, 
where the proton-nucleus and kaon-nucleus interactions are ignored. While 
the essential features of the high momentum transfer reaction (p, K + ) can be 
understood in this approach, the nuclear interactions may have some conse- 
quences. They produce both absorptive and dispersive effects [57]. However, 
for large incident energies considered in this calculation, the absorption effects 
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Fig. 12. Differential cross section for the 40 Ca(p, K + )\C&{Qd^i2) reaction for the 
incident proton energy of 1.5 GeV. various curves have the same meaning as in 
Fig. 11. 

are likely to be the most important; a qualitative estimate of this is provided 
in the following. 

We use the eikonal approximation to estimate the attenuation factor for a 
particle traveling through the nuclear medium. First, we define a refractive 
index of the nuclear medium as (see, e.g., [58]) 



r](r, E) = rj(E) 



p(r) 
Po : 



(39) 



where p{r)[p{\/b 2 + z 2 } is the nuclear density distribution, and 



V (E) = 



HE)' 



(40) 



Here k is the external wave number, and k is the wave number in the medium. 
The imaginary part of this refractive index is defined as r/o- In terms of 7/0, the 
attenuation factor can be written as [45] 
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Fig. 13. Angle integrated cross section for the 12 C(p, i^ + ) 1 lC(0p 3 / 2 ) reaction as a 
function of the incident proton beam energy. 



dh dz 



p(h,z)exp[-kr] (E)L(b)] 
J dh dz p(b, z) 



(41) 



where L{b) is the length of the path traveled by the particle in the medium, 
which is given by 



Lib) = '^dz (42) 
J Po 

If the nuclear density is approximated by a Gaussian function, 

p(r) = po exp (— r 2 /a 2 ), the integration in Eq. (41) can be done analytically. 

In this case the attenuation factor is given by 



The attenuation due to medium can be calculated if the value of 7]q(E) is 
known. This can be obtained from the imaginary part of the optical potential 

W as 
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Table 3 

Total attenuation (\S(E)\ 2 ) due to distortion effects in the initial and final channels. 



Reaction 


\S(E)\ 2 


4 He(p,^+)iHe(0s 1/2 ) 


0.28 


12 C(p,^+)^C(0p 3 /2) 


0.19 


40 Ca(p,K+) 4 ACa(0d 3/2 ) 


0.10 



Vo(E) = ¥ -W (E) (44) 
We use here the following high energy relations to obtain W 



_ 2 ka T p 



W (E)=h'^, (45) 

where ot is the total nucleon-nucleon or meson-nucleon cross section. In order 
to determine the total attenuation factor for the (p, K + ) reaction, the total 
attenuation due to both the proton and kaon distortions has been estimated 
by replacing the factor kr} in Eq. (43) by 



kr} -> k p r] (E p ) + k K r] (E K ) (46) 

In our estimation of the attenuation we have taken the average values of a? 
as 45 mb and 14 mb for pp [59] and K + p systems [60], respectively. The value 
of the parameter a is taken to be 1.61 fm, 2.37 fm and 3.52 fm for 4 He, 12 C 
and 40 Ca targets, respectively [61]. In Table 3, we show the values of the total 
attenuation \S(E)\ 2 for the three targets discussed in this paper. We see that 
overall effect of initial and final state interactions is to decrease the peak cross 
sections relative to the PW calculations by a factor of 3 - 10. We however, 
would like to emphasize that these estimates only provide an idea of the ab- 
sorptive effects. In more rigorous calculations, refractive distortion effects also 
play a role which could influence the shapes of the angular distributions as 
well. 



At the same time, one should also note that in the distorted wave treatment the 
continuum wave functions are no longer associated with sharp momenta but 
are states with a momentum distribution [see Eqs (30)-(31)]. This leads to a 
redistribution of the large momentum transfer differently from what is allowed 
in the plane wave approximation. It could shift the sensitivity of the model 
to even lower momenta leading to enhanced cross sections. The competition 
between this effect and the absorption effect would ultimately decide the role 
of distortions in these reactions. 
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Fig. 14. Angle integrated cross section for the m Ca.(p, K + ) 4 \Ca(0d 3 /2) reaction as 
a function of the incident proton beam energy. 



4 Summary and Conclusions 



In summary, we have made a study of the A(p, K + )^B reaction on 4 He, 12 C, 
and 40 Ca targets within a fully covariant general two-nucleon mechanism where 
in the initial collision of the incident proton with one of the target nucleons, 
AT* (1710), AT* (1650), and AT* (1720) baryonic resonances are excited which sub- 
sequently propagate and decay into the relevant channel. The initial nucleon- 
nucleon collisions are mediated by pion and also by rho and omega exchange 
mechanisms. Expressions for the reaction amplitudes are derived in both dis- 
torted wave (DW) and plane wave (PW) approximations. However, the numer- 
ical calculations are performed within the latter approximation which helps in 
understanding the essential features of the (p, K + ) reaction without requiring 
very lengthy and cumbersome computations which are necessarily involved in 
the DW method. Wave functions of baryonic bound states are obtained by 
solving the Dirac equation with appropriate scalar and vector potentials. 

In our model, the (p, K + ) reaction proceeds predominantly via excitation of 
the A/"* (1710) resonant state. For reactions on all the three target nuclei, the 
cross sections are dominated by those graphs in which the intermediate meson 
originates from the target nucleon and the projectile nucleon gets excited to 
the baryonic resonance state (target emission graph). The one-pion-exchange 



29 



processes make up most of the differential cross section at all angles, p and 
u exchange graphs individually are relatively more important at backward 
angles where the momentum transfer to the nucleus is comparatively larger. 
The effect of the heavier meson exchange is to reduce (although very weakly) 
the pion exchange only cross sections. However, the p and u exchange processes 
put together do not appreciably reduce the pion exchange only cross sections 
even at the back angles. 

The calculated cross sections are maximum for the hypernuclear state with the 
largest orbital angular momentum. This is the typical of the large momentum 
transfer reactions [57]. For heavier targets, the angular distributions for the 
favored transitions peak at angles larger than the 0° which in contrast to the 
results of most of the previous nonrelativistic calculations for this reaction. 
This reflects directly the nature of the Dirac spinors for the bound states 
which involve several maxima in the region of large momentum transfer. In 
case of the light target 4 He, however, the differential cross section still peaks 
near the zero degree as in this case the momentum transfers are in the region 
where the bound state spinors are smoothly decreasing with momentum. The 
energy dependence of the calculated total production cross section follows 
closely that of the pp — > pAK + reaction. 

The present results have been obtained using the PW approximation. Effects 
of initial and final state interactions have been estimated from an eikonal 
approximation for the nuclear elastic scattering interactions of the proton and 
kaon. The predicted cross sections are reduced by factors ranging between 
3 - 10. The reduction is less for the lighter targets. However, In the present 
treatment we considered the absorptive distortions effects only which influence 
the absolute magnitudes of the cross sections. In a more rigorous calculation, 
distortions may also affect the shapes of the angular distributions. 

Our results suggest that the study of the (p, K + ) reaction is attractive as 
it provides an alternative way to study the the spectroscopy of the A hy- 
pernuclear states. This reaction should be measurable at the COSY facility 
in the Forschungszentrum Jiilich. The characteristics of these cross sections 
predicted by us should be helpful in planning of such experiments. Such ex- 
periments may also lead to more quantitative calculations which should utilize 
the proper distorted waves. 
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A Expression for the amplitudes in distorted wave and plane wave 
approximat ions 



In this appendix we present the final form of the amplitude (M 2b ) in both full 
distorted wave as well as plane wave approximations. 



MgrW /2 ) = -^o(ff ) {9N l/2 N,){g Nl/2 KK + )- t 



1 



. 2m N ) WJ, i/2 J ' " > ^-1/2— > q 2 _ m 2 _ n ^( g ) 
1 



p%* - (m N * /2 -iY Nl/ j2) 2 _ 



£ EH* 

L v J v l K m K Mi 



x < L p 01/2fXi\J p fXi > ^2L p + lY eKmK (tt K ,tt' K ) 
c + ^E K )a WK Z^^^ (p A y 



~ <^W k Wk x - ip'K y WW> - ^)(3/2 + m) 
x v/(l/2 + *)(3/2 - riZ%£ffl^£Zl + i<PAr 

+ y K . zHx (pJ ) z7X;^ (p£> 

- ^^{(Pk, - ^ w )V(l/2-^)(3/2 + //i) 



x ^(1/2 + ^)(3/2 - ^)z^;(p^zTXr-l +1 (PAY 



where we have defined 



(A.l) 



(A.2) 



We have, c = -mA + m N *, and c = —1. Superscript (*) on a function 
represents the complex conjugate of the function. In Eq. (A.l) we have defined 



^(h, k 2 ) = 2Ag(k 1 )f(k 2 ) + Be(h)f(k 2 ) + Bgik^ih), 



(A.3) 
(A.4) 
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= 2Af(k 1 )g(k 2 ) + B£(ki)g(k 2 ) + Bf(h)e(k 2 ), 


(A.5) 




= fe A { k A) F L P J P {K, ki)ft K (k' K , k K ), 


(A.6) 




= ge' A (k A )F LpJp (k' i , h)fe K (k' K , k K ), 


(A.7) 




= fe A (kA)G UpJp (k , i , ki)fi K (k' K , k K ), 


(A.8) 




= ge' A (k A )G LpJp {k' i , ki)f eK {k' K , k K ). 


(A.9) 



In Eqs. (A. 4) and (A.5), A = and B — 1. There is one more similar term 
in Eq. (A.l) involving wave function ^{k A ) and i[i{k A ) i n place of fe(k A ) and 
ge(k\). In this term, coefficients Co and c' Q are replaced by c\ and c[ which are 
-I and 0, respectively. 

The form of the amplitude M 2b in the plane wave approximation is 



q 2 -ml- U n (q) 



J. 7T " \ / 



p%* - {m Nl/2 -iY N * / j2f 

+ z 7X l ^*^T 4 

+ <y +1 (PA) V(l/2-/i,)(3/2 + ^) T 2 

+ ^f +1 (PA)V( 1 /2-^)(3/2 + ^) T 5 
+ Z^f-\p A y ^(1/2 + ^(3/2- 



(A.10) 



where we have defined 



he 

Ti = (c + c E k ) f(k A ) - c'q— — Pip K sin 9 K f(k A ) 

hi + m N 

he 

+ {c l + c[E K ){ > {k A )c' 1 — PiPKCose K Z(k A ), (A.ll) 

hji + m N 

he 

T 2 = c' — — pip K sin 6 K f(k A ) 

hi + m N 

he 

+ c[—— PiPx sin9 K £(k A ), (A. 12) 

hi + m N 

T 3 = -T 2 , (A.13) 

he 

T 4 = -2i c' Q p K cos 6k g(k A ) + 2i -=— (-co + c' E K ) Pi g(k A ) (A.14) 

hi + m N 
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he 

- 2i c[ p K cos 9 K Z'(k A ) + 2i — (-ci + ci E K ) Pl £(k A ), (A. 15) 

T 5 = -ic' p K sin9 K g(k A ) ic[p K sin9 K £'(k A ), (A.16) 

T 5 = T 6 (A. 17) 

The coefficients c , Cq, Ci and c' x are the same as defined above. 

The amplitude for the graph 2(c) can be written in the analogous way. Those 
involving the excitation of other two resonances and vector meson exchanges 
have similar structure even though they may be a bit more complicated. 



B Continuum distorted wave functions in momentum space 



In this appendix we present a short derivation of Eqs. (30)-(31). 

The incident and outgoing particles are described by stationary continuum 
wave functions \f(r a ) = ip a (p a ,r a ) exp(—iE a t) where a — i represents the 
incident channel and K the outgoing one. ip a (Pa, r «) are given, in the partial 
wave representation, by 



4 +) (Pi, n) = £ iLp < l p m p - ^/MJpM P > Y^^ipi) 

J V L V M, 



M v 



[Ti 



pl/2J p \ 

P*(p K ,r K )= ]T r eK Ye K rnAPK)hAkK,a K )Y; KmK (r K ), (B.2) 
e K m K 



where various symbols have the same meaning as those described in the main 
text. Having chosen a momentum space formalism, we need to calculate the 
Fourier transforms 



*a( P ',Pa) = J ^- 4 e^ a (r a ), (B.3) 

depending on the unconstrained Fourier four momentum p' = (p' Q , p') and the 
asymptotic on-shell momenta p a = (E a ,p a ). The time integral in Eq. (B.3), 
can be trivially performed, leading to 



* a (p',p a ) = 5(p' -E a ) J |^e-^ r ^ a ( Pa ,r a ). (B.4) 
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Further integrations in Eq. (B.4) can be carried out by using the relation 



e-ip'-ra = 4 7r Y,i- L h(k'a a ) Y LM (p') Y£ M (r a ), (B.5) 

LM 

which lead to Eqs. (30)-(31). 
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